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Abstract Positronium is an ideal system for the research of the quantum elec- 
trodynamics (QED) in bound state. The hyperfine splitting (HFS) of positro- 
nium, Z\hfSi gives a good test of the bound state calculations and probes 
new physics beyond the Standard Model. A new method of QED calculations 
• has revealed the discrepancy by 15ppm (3.9tr) of Z\hfs between the QED 

prediction and the experimental average. There would be possibility of new 
physics or common systematic uncertainties in the previous all experiments. 
We describe a new experiment to reduce possible systematic uncertainties 
and will provide an independent check of the discrepancy. We are now taking 
data and the current result of Z\ H fs = 203.3951 ± 0.002 4(stat., 12 ppm) ± 
0.001 9(sys., 9.5 ppm) GHz has been obtained so far. A measurement with a 
. ' precision of O(ppm) is expected within a year. 
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1 Introduction 

Positronium (Ps), a bound state of an electron and a positron, is a purely 
leptonic system which allows for very sensitive tests of Quantum ElectroDy- 
namics (QED). The precise measurement of the hyperfine splitting between 
orthopositronium (o-Ps, l 3 ^) and parapositronium (p-Ps, l 1 5o) (Ps-HFS, 
^hfs) provides a good test of bound state QED. Ps-HFS is expected to be 
relatively large (for example compared to hydrogen HFS) due to a relatively 
large spin-spin interaction, and also due to the contribution from vacuum os- 
cillation (o-Ps — > 7* — > o-Ps). The contribution from vacuum oscillation is 
sensitive to new physics beyond the Standard Model. 



Experimental average Theory (20 00) 

Phys. Rev. Lett. 34, 246 (1975) 
Phys. Rev. A15, 241 (1977) 
Phys. Rev. A15 , 251 (1977 ) 
Phys. Rev. A27, 262 ( 1983) 
Phys. Rev. A30, 1331 (1984) 
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Fig. 1 Historical plot of zAjjfs- 

Fig. [1] shows the measured and theoretical values of Z\hfs ■ The combined 
value from the results of the independent and the most precise 2 experiments is 
A™£ s = 203.388 65(67) GHz (3.3 ppm) Q][l[3]. Recent developments in Non- 
Relativistic QED (NRQED) have added 0(a 3 In a) corrections to the theoreti- 
cal prediction which now stands at A$ FS = 203.391 69(41) GHz (2.0 ppm) [SEJ 
E]. The discrepancy of 3.04(79) MHz (15 ppm, 3.9cr) between A^ s and ^ h FS 
might either be due to the common systematic uncertainties in the previous 
experiments or to new physics beyond the Standard Model. 

There are two possible common systematic uncertainties in the previous 
experiments. One is the unthermalized o-Ps contribution which results in an 
underestimation of the material effect [7 . This effect has already been shown 
to be significant [8"|[9"|l 1 I I 1 l ^fiTlrTS"] in the o-Ps lifetime puzzle. Another is the 
uncertainty in the magnetic field uniformity which was cited as the most sig- 
nificant systematic error in the previous experiments. 

The energy levels of the ground state of Ps are shown as a function of static 
magnetic field in Fig. [5] Due to technical difficulties in directly stimulating 
^HFSi the previous experiments were indirect measurement by stimulating 
the transition of Zeeman splitting (Z\ m j X ) under static magnetic fields. The 
relationship between Z\hfs and ^Wx is approximately given by the Breit-Rabi 
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Fig. 2 Zeeman energy levels of Ps in its ground state. The arrows J, 4- means the spin up 
and down of electron, and the arrows -f|-, ij. means the spin up and down of positron. 

equation 

Zi mix ~ -A HFS (Vl + 4.T 2 - l) . (1) 

x is given as g' hbB/ (hAups), where g' = g (l — ^a 2 ) is the g factor for a 
positron (electron) in Ps [14,15,16,17]. hb is the Bohr magneton, B is the 
static magnetic field, and h is the Planck constant. Briefly, Z\ m i x is measured 
by applying microwaves under the static magnetic field and then Z\hfs is 
obtained from Eq. ([1]). The experimental signal is the change of 27 and 37 
annihilation rate caused by Zeeman transition. The experimental resonance 
line shape is obtained as explained in the ref. [7]. 

2 Experimental Setup 

The photograph of our new experimental setup is shown in Fig. [3] and the 
schematic diagram of the setup is shown in Fig.UJ The characteristic features of 
our experiment are the large bore superconducting magnet, /^-tagging system, 
and high performance 7-ray detectors. Details are discussed in the following 
sections. 



2.1 Large bore superconducting magnet 

A large bore superconducting magnet is used to produce the static magnetic 
field B ~ 0.866 T, which induces the Zeeman splitting. The bore diameter 
of the magnet is 800 mm, and its length is 2 m. The magnet is operated in 
persistent current mode, making the stability of the magnetic field better than 
±lppm. With the compensation coils, the uniformity of the magnetic field is 
0.9 ppm (RMS) in the large volume of cylinder 40 mm in diameter and 100 mm 
long, where Ps are formed. 
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Fig. 3 Photograph of the setup of our new experiment. The large bore superconducting 
magnet and the microwave waveguide are shown. Microwaves are guided through the waveg- 
uide into the cavity. The RF cavity, /3-tagging system, and the 7-ray detectors are located 
at the center of the magnet. 
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Fig. 4 Schematic diagram of the setup (top view in magnet). 



2.2 /3-tagging system 

The positron source is 14 [id (500 kBq) of 22 Na. A plastic scintillator 10 mm 
in diameter and 0.2 mm thick is used to tag positrons emitted from the 22 Na. 
The scintillation light is detected by fine mesh PMTs and provides a start 
signal which corresponds to the time of Ps formation. The positron stop in 
the microwave cavity, filled with 1-C4H10 gas, forming Ps. 

Ps decays into photons that are detected with LaBr3 (Ce) scintillators. The 
timing difference of positron emission and 7-detection is used to improve the 
accuracy of the measurement of Z\hfs a $ follows: 
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1. Imposing a time cut means that we can select well thermalized Ps, reducing 
the unthermalized o-Ps contribution. 

2. A time cut also allows us to remove the annihilation events (prompt peak 
in Fig. Eta)). The S/N is improved significantly by about a factor of 20. 

2.3 High performance 7-ray detectors 

Six 7-ray detectors are located around the microwave cavity to detect the 
511 keV annihilation 7-rays. LaBr3 scintillators, 1.5 inches in diameter and 

2 inches long are used. LaBr3 scintillators have good energy resolution (4% 
FWHM at 511 keV) and timing resolution (0.2ns FWHM at 511 keV), and 
have a short decay constant (16 ns). 

3 Analysis 

The experiment was performed from July 2010 to 11th March 2011. Currently 
the run was stopped by the earthquake occurred on 11th March 2011, but will 
be restarted soon. In the overall periods, the trigger rate was about 2.4 kHz 
and the data acquisition rate was about 1.2 kHz. The data acquisition was 
performed using NIM and CAMAC systems. 

Fig. [5ja) shows one of the measured timing spectra. The prompt peak con- 
tains the two processes, annihilation and |— } decay. Decay curve of |+) and 
m z = ±1 states contribute dominantly after 30 ns, then the constant acci- 
dental spectrum. A timing window of 35-155 ns is applied to select transition 
events. Fig. [3£b) shows the energy spectra after the timing cut. It is obtained 
by subtracting the accidental contribution using the timing window of 950- 
1350 ns. 

The difference between RF-ON spectrum and RF-OFF spectrum is the 
Zeeman transition amount essentially. The detection efficiency which is neces- 
sary to obtain the transition amount is estimated by Monte Carlo simulation. 
The 27 and 37 decay rate were measured at various magnetic field strengths 
with a fixed RF frequency as shown in Fig. [5] Fig. [5] shows one of the reso- 
nance lines we obtained. The data are fitted by theoretical resonance line. To 
estimate the material effect, the resonance curves have been measured at four 
gas densities. The fitting results are summarized in Table [T] 

The gas density dependence of Z\hfs is shown in Fig. [7] Currently Z\hfs 
at vacuum is obtained by extrapolating the data linearly. The fitting results in 
zIhfs = 203.395 1 ± 0.002 9 GHz (14 ppm) and the gas density dependence is 
—248 ±21 ppm/amagat. The error includes statistical uncertainty, uncertainty 
from statistics of Monte Carlo simulation, uncertainty of RF power, and un- 
certainty of Ql value of the RF cavity. Systematic errors of the current result 
are summarized in Table [5] 

The current result is 



zAhfs = 203.395 1 ± 0.002 4(stat., 12 ppm) ± 0.001 9(sys., 9.5 ppm) GHz , (2) 
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Fig. 5 Timing and energy spectra. The black lines are RF-ON, and the gray lines are RF- 
OFF. (a) Decay curves of Ps. The decay rate of Ps increases with RF because of the Zeeman 
transition, (b) Energy spectra. 
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Fig. 6 Resonance line at 0.895 15(20) amagat gas density. The circles and error bars are the 
data, and the solid line is the best fit result. The lower figure shows the ratio of measured 
data over best-fit function. The error bars include errors from statistics of data, statistics of 
Monte Carlo simulation, uncertainty of RF power, and uncertainty of Ql value of the RF 
cavity. 



Table 1 Current fitting results of the resonance lines. These uncertainties include errors 
from statistics of data, statistics of Monte Carlo simulation, uncertainty of RF power, and 
uncertainty of Ql value of the RF cavity. 



Gas density (amagat) 


zAhfs (GHz) 


Relative error (ppm) 


X 2 /ndf (p) 


0.895 15(20) 


203.350 6(20) 


9.8 


4.1/5 (0.54) 


0.673 99(19) 


203.359 3(25) 


12 


6.7/4 (0.15) 


0.23491(11) 


203.384 6(24) 


12 


3.7/3 (0.30) 


0.168 434(21) 


203.383 3(51) 


25 


2.7/4 (0.61) 
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Fig. 7 Gas density dependence of Z\hfs ■ The circles and error bars are the data, the solid 
line is the best-fit with linear function, the blue band is la uncertainty, the green band is 
experimental average at vacuum, and the red band is 0(cr 3 In a) QED prediction. 



Phys. Rev. Lett. 34, 246 (1975; 
Phys. Rev. A15, 241 (1977) 
Phys. Rev. A15, 251 (1977) 
Phys. Rev. A27, 262 (1983) 
Phys. Rev. A30, 1331 (1984) 
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Fig. 8 Our current result in historical plot of Z\hfs ■ 



Table 2 Summary of systematic errors of our current result. 



Source Errors in zIhfs (pP m ) 



Magnetic Field: 



Non-uniformity 


1.8 


Offset and reproducibility 


1.0 


NMR measurement 


1.0 


Detection Efficiency: 




Estimation using simulation 


7.0 


Material Effect: 




Ps thermalization 


3.0 


RF System: 




RF power 


2.9 


Ql value of RF cavity 


4.3 


RF frequency 


1.0 


Quadrature sum 


9.5 



which is shown in Fig. [5] compared with previous experiments and theoretical 
calculation. 
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4 Prospects 

1. Material Effect Currently we assumed that Z\hfs depends on gas den- 
sity linearly. If the unthermalized Ps contribution is large, the depen- 
dence becomes non-linear. According to the previous measurement of the 
momentum-transfer cross section [TB] , unthermalized effect is estimated to 
be less than 3ppm with 1-C4H10 gas. We are now precisely measuring the 
Ps thermalization using the same technique as rcf. 8 , 9 [ [TOlfTTlfT^lfTB] (pick- 
off), whose technique is different from that of ref. [TB] (Doppler-broadening 
spectroscopy). 

2. RF System We will carefully control the experimental environment, espe- 
cially the temperature, which will reduce the uncertainty. 

3. Detection efficiency Currently it is estimated by Monte Carlo simulation. 
It will be carefully studied and will be estimated by real data, which also 
will reduce the uncertainty. 

4. Statistics I2ppm has been obtained. By measuring more efficient points of 
gas density and magnetic field strengths, a measurement with a precision 
of O(ppm) is expected within a year. 

5 Conclusion 

A new experiment to measure the Ps-HFS which reduces possible common 
uncertainties in previous experiments has been constructed. The current result 
of Z\ H fs = 203.395 1 ± 0.002 4(stat., 12 ppm) ± 0.001 9(sys., 9.5 ppm) GHz has 
been obtained so far. A new result with an accuracy of O(ppm) will be obtained 
within a year which will be an independent check of the discrepancy. 
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